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We measure the resonance line shape of atomic vapor layers with nanoscale thickness confined 
between two sapphire windows. The measurement is performed by scanning a probe laser through 
resonance and collecting the scattered light. The line shape is dominated by the effects of Dicke 
narrowing, self-broadening, and atom-surface interactions. By fitting the measured line shape to a 
simple model we discuss the possibility to extract information about the atom-surface interaction. 


I. INTRODUCTION 

Atomic alkali vapor cells are widely used in applica¬ 
tions ranging from magnetometry [T], sources of quan¬ 
tum light [2], electric field imaging [3], atomic clocks [4], 
nanophotonics [5], optical isolators [6] and quantum 
memory [7]. Modern vapor cell technology has been 
trending towards integrating miniaturized cells into the 
aforementioned applications, with millimeter iHii] , mi¬ 
cron m and even nanometer scale cells m popular ar¬ 
eas of investigation. An accurate method of analyzing 
spectra from miniature cells is vital to aid the develop¬ 
ment of such technology. 

If the vapor is confined spatially with a dimension less 
than the transition wavelength fascinating new physical 
effects become accessible such as the narrowing of spec¬ 
tral lines la, extreme dispersion resulting in large nega¬ 
tive group indices and superluminal propagation [Tgns], 
repulsive van der Waals interactions m, the coopera¬ 
tive Lamb shift m, and perhaps a medium where the 
Kramers-Kronig relations are violated m- 

The tight confinement of the atoms inside nanocells 
opens up opportunities to study the interaction of atoms 
with a nearby surface, explored using spectroscopy on 
both low-lying m and higher-lying excited states [21] or 
EIT spectroscopy of highly-excited Rydberg states [22] . 
This can be expanded to investigate the temperature de¬ 
pendence of the coefficients describing the strength of the 
atom-surface (AS) interaction [23] and cases where the 
usually attractive AS interaction becomes repulsive due 
to surface resonances m- Alternative methods to spec¬ 
troscopy for probing the AS interaction have also been 
used, such as scattering or deflection of a beam close to a 
surface [24ll29] . atomic beam diffraction [301131] and re¬ 
flection of an ultracold atom cloud from an atomic mir¬ 
ror [32l [33]. Such experiments use detection methods 
that take place some time after the interaction has oc¬ 
curred. 

In this paper we employ spectroscopic measurements 
that allow us to infer the effects of the AS interaction 
at the time of absorption. Note that in this case, both 
ground and excited states contribute to the measurement 
of the AS interaction. In what follows we present de¬ 
tailed spectra for a range of cell lengths and tempera¬ 


tures, expanding upon a previous publication [20], giv¬ 
ing further detail into the error analysis and fitting pro¬ 
cedures, present a simulation justifying the model used 
to model velocity selective effects used and offer further 
analysis on the results found. We describe an experi¬ 
ment that takes many absorption spectra over a range of 
temperatures and cell lengths, and uses fitting and error 
analysis to precisely measure the general form of the AS 
interaction within the near field. We first detail the meth¬ 
ods of data acquisition for absorption spectra, where we 
utilize single photon counting modules (SPCMs) to ac¬ 
quire absorption spectra over long integration times. The 
ease of acquisition allows us to take many spectra for a 
range of cell lengths and temperatures. We then discuss 
how experimental error is accounted for and utilized to 
increase the precision of our measurements. 

The theory section details the form of the AS interac¬ 
tion in the near field (van der Waals) regime [34], dis¬ 
cussing possible effects that may alter the final form, 
such as a transition between the near field and retarded 
(Casimir Polder) [35] regimes; the effect of multiple re¬ 
flections and the effect of temperature on the interac¬ 
tion. We then give an overview of the fitting procedure 
performed on the spectra taken. The model we have de¬ 
veloped accounts for self broadening [36] , Dicke narrow¬ 
ing effects [H] and the AS interaction [34]. The results 
section presents fitting results for spectra taken on the 
Rb D2 line and the Cs D1 line, where we find that to 
describe spectra in the length range investigated herein, 
the van der Waals description of the AS interaction is 
optimal. We finally discuss the implication and accuracy 
of our results including a discussion of the limitations of 
the technique and perspectives for further work. 


II. EXPERIMENT 
A. Data Acquisition 

To probe the atom surface (AS) interaction we use 
nanocells (NCs), produced at the National Academy of 
Science in Armenia m- The cells are made from c-axis 
aligned sapphire with an interior surface roughness of 
3 nm and are described in more detail in previous pub- 
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FIG. 1. Experimental schematic used to measure the atomic 
line shape. The laser is scanned through resonance and we de¬ 
tect the scattered light. Part of the beam is split off by a half¬ 
wave plate (A/2) and polarizing beam splitter cube (PBSC) 
to the reference optics for frequency calibration. The beam is 
then hbered to the nanocell (NC) section of the experiment, 
where lens LI focuses the beam to a 1/e^ radius of 20 /xm 
inside the NC. The transmitted beam is re-collimated by lens 
L2 and a transmission spectra is collected on photo-diode PD. 
Off-axis scattered photons are collimated by lens L3, with am¬ 
bient light filtered out by an infrared (IR) hlter, and sent into 
a single photon counting module (SPCM) via an hber colli¬ 
mator FCl. The signal is processed by the oscilloscope to 
generate a histogram of the photon arrival times. 


lications [mm EH Eg. Each cell has a wedged profile, 
with a thickness range of 2 /xm to 30 nm. In order to 
exploit variable vapour column thickness in the range of 
30-2000 nm, the cell is vertically wedged by placing a 
2 /xm thick platinum spacer strip between the windows 
at the bottom side prior to gluing. 

The NC thickness, L, is measured using back reflec¬ 
tions from the surfaces inside the cell. Two back re¬ 
flections Rl, from the front surface of the cell, and R2, 
a combination of reflections from both interior surfaces, 
are measured. Etalon effects inside the cell mean the 
two beams that compose R2 interfere constructively or 
destructively depending on the cell length, resulting in a 
varying intensity of R2. The thickness is calculated by 
taking the ratio of intensities Rl and R2 and applying a 
standard Fabry-Perot treatment [40]. 

Our method measures absorption or more precisely ex¬ 
tinction in the forward scattering direction by detecting 
the off-axis scattered light resulting in absorption spectra 
that we can fit to measure the AS interaction. We use 
this method because a measurable transmission profile 
for extremely short cell lengths can only be obtained by 
heating to high temperatures (exceeding 300°C), where 
collisional broadening smears out potential atom surface 
measurements [36|. For shorter cell lengths (L < 100 nm) 
and lower temperatures (T < 200°C), absorption is less 
than 0.5% of the incident light, making the measurement 
of transmission impractical. Using off-axis detection of 
scattered light with single-photon counters, extremely 
high sensitivity spectra can be taken from shorter cell 
lengths at lower temperatures. This increases our exper- 
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FIG. 2. Off-axis scattered photons measured using an exci¬ 
tation laser scanned across the Rb D2 resonance. The figure 
shows a comparison of raw data (gray areas, panels (a) and 
(b)) to processed data (black points, panels (c) and (d)) on 
the Rb D2 line at a temperature of 125°G for (a,c) 60 nm 
and (b,d) 250 nm cell lengths. Background counts are shown 
in gray and the atomic signal highlighted in red and blue 
for lengths of 60 and 250 nm respectively. Panel (a): The 
atomic signal, highlighted in red, is small in comparison to 
the background count because the interrogation volume is so 
small that there are relatively few atoms contributing to the 
signal. Therefore, the error is dominated by the error of the 
background count. Panel (b): The atomic signal is the main 
component of the raw signal because far more atoms are con¬ 
tributing to the signal, resulting in much smaller relative error 
compared to panel (a). Panel (c): The ht for L = 60 nm (red 
line) has a ‘goodness of ht’ parameter, T, of 0.1. The small 
value of T is caused by the large background counts involved 
for such small cell lengths, illustrated in panel (a). Panel (b): 
the ht for L = 250 nm (blue line) has a reduced chi squared 
T [33 of 1.7, indicating an excellent ht. The background pho¬ 
ton count is much smaller than for L = 60 nm, resulting in 
lower relative error and an T closer to 1. 


imental resolution as spectra with all hyperhne states re¬ 
solved can be taken using single-beam spectroscopy. The 
narrowing of the lines is caused by Dicke narrowing [T4| 
effects inside the cell- a motional effect that will be out¬ 
lined in section [ml Absorption spectra were taken on 
the Rb D2 line; used for its large AS interaction on the 
ground state transition, and the Cs D1 line; used for its 
large 9 GHz ground state splitting and 1.2 GHz excited 
state splitting, much larger than the Doppler width. All 
spectra were taken within the weak-probe limit m- 

The experimental layout is shown in Fig. Some of 
the light is sent though reference optics - a 7.5 cm ref¬ 
erence vapour cell and a Fabry-Perot etalon, which are 
collected concurrently with the scattered light spectra to 
calibrate the laser frequency. The light is then sent down 
an optical fiber to the NC section of the experiment. The 
light is focussed in the centre of the NC by lens LI to a 
spot size with a 1/e^ radius of 20 /xm, then is recollimated 
by lens L2. Transmission data are recorded on a photodi- 
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ode (PD). Off-axis scattered light is collected by lens L3, 
with background light and thermal photons from the cell 
heater hltered out by an infrared hlter. The signal is sent 
to a single-photon counting module (SPCM) via a hber 
collimator, FC3. The counts are processed by a LeCroy 
Waverunner 610Zi oscilloscope into a histogram of arrival 
times corresponding to the frequency of the laser scan, 
according to the method outlined in [38] . 

The laser scan time is 50 ms, and individual bins in the 
histogram are 5 /is long, corresponding to a bin width of 
2 MHz. Hence the detection method is sensitive to the 
frequency at which a photon was absorbed, meaning AS 
shifts measured map directly onto a particular AS separa¬ 
tion for that frequency. Our method is insensitive to any 
motion between the time of absorption and fluorescence 
events. Single scans return only a few photon counts, so 
long integration times (15 minutes to several hours) are 
used to build up statistics of arrival times, allowing for 
sensitive detection of spectra. 

Figs, [^a) and (b) show examples of raw spectra taken 
for T = 125°C. The data are taken at 2 different cell 
lengths, 60 nm (red area) and 250 nm (blue area). The 
raw spectra have a very large background count, high¬ 
lighted in grey. Signal-to-background ratio varies hugely 
between different cell lengths, only the atomic signal is 
highlighted (colored areas) in Figs, [^a) and (b). The 
effect of these background counts on experimental errors 
will be discussed in the next subsection. Before analy¬ 
sis, the background signal is subtracted from the atomic 
signal. The resulting processed data is shown as black 
points in Figs, [^c) and (d). We ht the spectra, using a 
procedure outlined in section [111} Examples of ht results 
are shown on Figs, [^c) {L = 60 nm, red line) and (d) 
(L = 250 nm, blue line). 


B. Error analysis 

The errors in the atomic signal, A{uj)^ are calculated 
by taking the frequency-dependent raw spectrum R{uj) 
and subtracting the background count B. 

A{uj) = R{uj) - B. (1) 

For data acquisition, the errors in the atomic signal (sa) 
are Poissionian counting statistics; each data point in 
the raw counts R{uj)^ with a number of counts Account, 
has an associated error of VAccount [32]. The errors can 
be calculated as: 


— \J s\ + s|, (2) 

with Sr and sr as the errors in the raw and background 
signals, respectively. The relative error bars are very 
small for spectra similar to Fig. because A{uj) is large 
and B is small. However for spectra with low signal-to- 
background ratio like that in Fig. ^a), i.e. B ^ A{uj)^ 
the error of the raw signal is roughly equal to the error 


of the background sign al, s r ^ sr. Hence, sa can be ap¬ 
proximated as sa ~ \/2s|. This results in large relative 
error for shorter cell lengths when compared to the er¬ 
ror bars for longer cells with better signal-to-background 
ratios. 

We ht the spectra using the Marquardt-Levenburg 
method [37] to minimise the sum of the squares of the 
difference between theory and experiment, normalized by 
the error bar. For errors in htted parameters for an in¬ 
dividual spectrum, we use the standard deviations of the 
mean calculated by the htting routine. 

For the analysis performed on Rb spectra, we take an 
average of the htted parameters weighted by the reduced 
X-squared of the ht. The reduced y-squared is a good¬ 
ness of ht measurement, represented here with T to avoid 
confusion with the vapor susceptibility. T is calculated 
as 


T 


1 [A{uj) - T{u)]‘^ 

7 / / ^ 


(3) 


where A{uj) is the observed frequency-dependent atomic 
signal, T{u) is the theoretical signal and u is the degrees 
of freedom, calculated as = Appoints ~ A^param, where 
Appoints is the length of the data set and A^param is the 
number of htted parameters. 

Values for T of our htted spectra vary between 0.1 to 
over 1000. Generally, a T close to 1 indicates a good 
ht, and any value much lower than 1 indicates overesti¬ 
mated errors [37]. However, this is not the case in our 
experiment; the smaller values of T are caused by the 
huge background counts for shorter cell lengths, illus¬ 
trated in Figs. I^a) and (b). The small values of T can 
be accounted for by considering the large relative error 
margins that arise when the atomic signal is small result¬ 
ing in a higher proportion of the htted spectra within the 
error margins of the atomic signal. Consequently, the re¬ 
sulting T is much less than 1, but still a valid measure 
of goodness of ht. 


III. THEORY AND FITTING 
A. The Atom-Surface Interaction 

The atom-surface (AS) interaction is the attraction be¬ 
tween a dipole and its rehection in a surface. An in¬ 
duced dipole in the atomic medium induces a shift of 
charge in a nearby conducting or dielectric surface [34]. 
This charge redistribution results in a virtual dipole im¬ 
age [35] , which has an attractive interaction with the orig¬ 
inal atomic dipole. This causes a red-shift of the atomic 
transitions, generally described using —Cc^/r^, where a 
is an exponent that varies between 3 and 4, depending 
on the interaction regime. 

The Casimir-Polder (CP) interaction describes this 
process at AS separations r > A/27r, where A is the 
transition wavelength. The interaction between the real 
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FIG. 3. (a) The position-dependent red shift from the vdW 
1/r^ AS potential. Inset: The effect of multiple reflections 
contributing to the vdW AS interaction. Plotted are shifts 
caused by accounting for a single (blue line); two (purple dot¬ 
ted line); three (red dashed) and four reflections (green dot 
dashed line). These differences are on the order of 100 kHz, 
smaller than the range of shifts our experiment is sensitive to. 
Panel (b) shows the length-dependent variation of AS poten¬ 
tial (purple line) taken from equation (45) in |42] , highlighting 
the transition between the short range van der Waals and long 
range Casimir-Polder regimes. At short lengths, there is good 
agreement with a 1/r^ vdW potential (red dashed line). At 
longer lengths there is reasonable agreement with a 1/r^ CP 
potential (blue dot dashed line). U^dw is scaled such that 
Fvdw = 1 at f = A/27r, and f is the AS distance in units of 
X/27T. Most of our data come from the region highlighted in 
gray (length range of 10-40 nm). In panel (c), the variation 
of the exponent a of the atom-surface attraction in the form 
1/r^ between the vdW and CP regime is shown, again with 
the length range the experiment is sensitive to highlighted in 
gray. 


and image dipole is retarded by the transit time of 
photons [34], and the potential takes the form Ucp = 
with (^4 as a coupling coefficient determining 
the strength of the interaction. 

When the AS separation reaches the regime r < A/27r, 
the induced shifts follow the van der Waals description 
t^vdw = —C's/r^, where C 3 is the coupling coefficient 
of the van der Waals interaction. Fig. |^a) shows the 
induced red shift caused by t/vdw- 

Cd is dependent on the atomic polarizability m, and 


hence varies depending on the electronic transition the 
shift is acting upon. When measuring the AS shift on 
an optical transition, the effective coupling coefficient is 
the difference between the ground and excited states. 
For this reason, we distinguish between coupling coef¬ 
ficients using a superscript denoting the transition, e.g. 

for the ^^Rb bSi /2 Fg = 2,3 ground states and 
^5 S'i/ 2^5P3/2 measurements taken on the ^^Rb 55 'i /2 
Fg = 2,3 to ^^Rb 5 P 3/2 Fe = 1,2,3,4 transition. 

The smooth transition between the vdW (near field) 
and CP regime [43] has been the subject of many stud¬ 
ies [35] [42j [44] , with estimates for the onset of the transi¬ 
tion ranging from lengths of A/27r [42] to 0.03A [45]. To 
identify the length regime our work covers, we inspect 
the full form of the AS interaction including both CP 
and vdW forms [42] : 


Uas = - 




j,4l 

jFJa. 


-1 


(4) 


where M is the atomic mass and are length parame¬ 
ters related to the strength of the vdW/CP potential and 
are taken from [42] . Values taken are for ground state Rb 
(55 'i/ 2), unlike the transitions we investigate experimen¬ 
tally, which probe the difference in interaction coefficients 
between the ground (55'i/2) and excited ( 5 P 3 / 2 ) states. 
Calculation of /33^4 for the excited state transitions are 
beyond the scope of this study, and we use the informa¬ 
tion from Fig. [^b) as an indicator for which interaction 
regime we expect our spectra to lie in. We expect that the 
ground-state case is a worse case scenario as we are fur¬ 
ther into the Cs region for the excited state as the tran¬ 
sition wavelengths to nearby states are longer. Fig. I^b) 
shows the resulting potential for the ground state (pur¬ 
ple line). The length scale f is scaled in units of A/27r. 
Pvdw is normalized such that Pvdw = 1 at f = 1. Both 
the CP and vdW interactions are plotted for comparison 
(blue dot-dashed and red dotted lines respectively). The 
length range that most of our data are produced for is 
shown as a gray area on the plot. In Fig.j^c) we plot the 
length dependence of the exponent a of the AS interac¬ 
tion. Figs, [^b) and (c) demonstrate that our experiment 
is sensitive to the shorter range van der Waals interac¬ 
tion regime. The average exponent in the length range 
P = 10 — 40 nm is a = 3.16. 

The impact of the AS interaction on the final reso¬ 
nance line shape is usually an overall red-shift of the line 
shape, accompanied by a long asymmetric red tail due to 
atoms at different positions inside the cell experiencing a 
different shift. However, if surface polariton resonances 
coincide with the atomic resonance, the AS interaction 
becomes repulsive, instead causing a blue shift of the 
atomic resonance m- 

The values of C 3 measured in this experiment are for 
the Rb 55 'i/2 ^ 5 P 3/2 transition. The Cs coefficient for 
the transition when the atoms are interacting with a per¬ 
fectly conducting surface is = 4.1 ±0.2 kHz 

/im^ m, calculated using the methods outlined in m- 
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When the surface is not a perfect conductor, C 3 is mul¬ 
tiplied by a reflection coefficient R [48] . Sapphire has no 
significant variation in the frequency-dependent permit¬ 
tivity e(co’) in the relevant frequency range, and hence we 
use a value of esapph(<^) = 3.24 for the Cs D 1 and Rb 
D2 line, hence the reflection coefficient can be expressed 
as R = (esapph(^) - l)/(esapph(^) + !)• This results in a 

scaling factor of 0.53 for values of measured 

in this experiment. 

The coupling coefficient C 3 also varies with tempera¬ 
ture [23l [43] . However at short ranges the thermal de¬ 
pendence of Cs is dominated by surface excitations [23] . 
This is not the case for the transitions investigated herein, 
as they do not concur with the frequencies of any sur¬ 
face plasmon excitations of sapphire [48]. There is still 
some temperature variation of the Cs coefficient, this is 
on the order of 1 % over the temperature ranges investi¬ 
gated herein and is hence negligible in the scope of this 
experiment [46] . 

The AS interaction is also affected by multiple reflec¬ 
tions in the surface, where the induced dipole interacts 
with reflections at integer multiples of the AS separation 
r. The effect of such reflections are shown in the inset to 
Fig. 1^. The potential from a single reflection (blue line) 
differs only slightly from the effect of two (purple dot¬ 
ted line), three (red dashed) and four (green dot-dashed) 
reflections. These differences are on the order of a few 
hundred kHz, smaller than detectable in the experiment. 


B. Line shapes for absorption spectra in 
nanometric cells 


We analyze the raw data using a model developed to 
describe the susceptibility of an alkali vapor. The fitting 
function calculates the detuning-dependent susceptibil¬ 
ity x(A) of the vapor, then converts it to an absorp¬ 
tion spectrum, using a process similar to that outlined 
in [49] . The elements included in the model are outlined 
in Fig. [^ Each step is illustrated with an example line 
shape with the most relevant equation displayed beneath 
it. For comparison, the state of the line shape after (be¬ 
fore) each step is shown in red (gray). 

The first step is calculating an initial spectrum contain¬ 
ing the natural line-width, Fq, for each hyperfine transi¬ 
tion, shown in panel (a) on Fig. Each transition is an 
independent 2 level system. We then sum over all tran¬ 
sitions. The single 2-level susceptibility is given by m- 


x(A) = 


iK 

ro/2-iA’ 


(5) 


where A is the detuning and is a proportionality con¬ 
stant that includes the dipole matrix element and tran¬ 
sition strengths. Atomic data including the natural line- 
width, dipole matrix element, transition strengths and 
level splittings are taken from [49] . 

Panel (b) of Eig.j^shows the effects of collisional broad¬ 
ening, a density-dependent broadening caused by inter¬ 


atomic collisions. When an atom in the excited state 
undergoes a collision with another atom, there is a prob¬ 
ability of the excited atom returning to the ground state. 
This results a reduction of the lifetime leading to a homo¬ 
geneous broadening of the spectral lines that is linearly 
dependent on density M [50]. We model the collisional 
broadening Fcoi as Fcoi = where /3 is a proportion¬ 
ality constant, using values from [36]. The density is cal¬ 
culated from the measured temperature using the vapor 
pressure formulas in [49]. We assume we are within the 
collision regime where the binary approximation is valid. 
This is accounted for in the model by including it in the 
calculation of the susceptibility for single transitions, re- 
placing To in Eq.|5]with Ttot, where Ttot = ro+Fcoi. The 
line shape is then summed over all transitions to generate 
a spectrum describing the natural line shape broadened 
by inter-atomic collisions. 

Next, we account for velocity-dependent effects, illus¬ 
trated in panel (c) of Eig. In longer cells, the un¬ 
restricted motion of atoms in all directions leads to a 
Doppler broadening of the spectral lines. However, the 
tight confinement of the NC restricts the motion of atoms 
in the propagation direction, altering velocity-dependent 
shifts seen in the spectra. Atoms with a low velocity 
vector in the direction of laser beam propagation have a 
longer interaction time and therefore dominate the sig¬ 
nal. 

The velocity selected line shape can be understood by 
exploring the effect of the cell geometry on the atom- 
light interaction, illustrated schematically in Eig. ??(a). 
The NC is much shorter than the mean free path of 
the atoms, hence atoms traveling perpendicular to the 
windows with a high velocity in the direction of beam 
propagation, Vz^ collide with the opposite wall frequently, 
suppressing their contribution to the signal m- Atoms 
traveling at a small angle 0 to the wall with a low Vz 
have a much longer interaction time, and a much smaller 
Doppler shift. These ‘slow’ atoms undergo more ab¬ 
sorption relative to the ‘fast’ atoms, resulting in en¬ 
hanced sub-Doppler features on top of a small Doppler- 
broadened background in NC spectra. 

The extreme narrowing is illustrated in Eig. ??(b), 
where the narrowed line shape (purple line) is compared 
to the Doppler broadened absorption spectrum expected 
for longer cell lengths (gray area). At low temperatures 
(< 100 °C), this narrowing can be so extreme that all in¬ 
dividual hyperfine states can be resolved on the Rb D 2 
line (see Eig. 2(a) in [20]). 

The model accounts for such velocity effects using a 
phenomenological bimodal velocity distribution for the 
density of atoms, first described in [38] It describes the 
number of atoms with velocity v, as: 

Nyc{v) = C{aGf^st{v) + (1 - a)Gsiow(^)), (6) 

where C = {au^/7^ -h (1 — a)fu^/7^)~^ is a normalization 
constant, with u = y^ 2 /cbT/M, the rms velocity of the 
atoms, and a is a coefficient describing the strength of 
the contribution of traditional Maxwellian velocities to 
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FIG. 4. Illustration of the various contributing factors to the measured line-shapes. From left to right, the line-shape after each 
effect (red line) is added to the previous line shape (shown as grey area). Each step is accompanied by the most relevant equation 
used to describe the changes in the line shape. The spectra are generated for a cell length of 60 nm and temperature of 200°C. 
Panel (a) shows the natural line shape at zero velocity, — 0)- panel (b), effects caused by density (N) dependent self 

(collisional) broadening are added to the widths of each transition. This results in a large Lorentzian broadening of the line. 
Panel (c) shows velocity-dependent effects on the line shape. The Doppler-broadened line shape x(A,Uz) is convolved with 
a bimodal velocity distribution g(v) that represents the stronger contribution of ‘slow’ atoms to the line shape. The effect is 
a small additional width. In panel (d) the AS interaction is accounted for by adding a shift modeled using the AS potential 
fAdw for the ground to excited state acting over the entire cell length. The change on the line shape is dramatic for short cell 
lengths, shifting the peak and causing a long red tail. 


(a) 




FIG. 5. (a) Atomic trajectories in the nanocell that cause 
Dicke narrowing. The cell is shorter than the mean free path 
of the atoms, meaning that atoms traveling perpendicular to 
the cell walls, with a high velocity in the z direction (vz) 
do not have enough time to contribute to the signal, leaving 
atoms traveling at a small angle to the wall (0) with a much 
smaller Vz as the main component of transmission and absorp¬ 
tion spectra, (b) A comparison of Dicke-narrowed absorption 
spectra (purple line) to the Doppler-broadened spectra (gray 
area) expected in a 7 cm cell, highlighting the dramatic nar¬ 
rowing caused by velocity-selective effects inside the cell. 


the signal, and is one of the floating parameters in the 
fit. Nyc{v) is split into 2 components, Gfast('^)- the usual 
Maxwell distribution for describing the Doppler broad¬ 
ened profile, weighted by fitting coefficient a: 

Gfast(v) = exp(-v^/u^). (7) 

Gsiow('^) describes a narrowed Gaussian profile, repre¬ 
senting the atoms traveling at small angles to the walls, 
that make the larger contribution to the signal: 


( 8 ) 



Distance from cell walls (nm) Velocity (km/s) 


FIG. 6. Monte-Garlo simulations of atomic trajectories to test 
the appropriateness of the Dicke narrowing model used in our 
fitting. Simulations are performed with no AS force present 
(blue lines) and including the effect of the AS force (dashed 
red lines). These are compared to the bimodal distribution 
used in our fitting procedure (thick black lines). In panel (a) 
the weighting of contributions to the line shape for atomic 
positions inside the cell is shown. The gray dashed lines mark 
10 nm from the cell walls, where the AS force starts to have a 
considerable effect on the atoms. Panel (b) shows the velocity 
distribution of atoms. The distribution used in our model 
(black line) is in reasonable agreement with the simulations 
for atoms unaffected by the AS force (blue line) and atoms 
which have been accelerated towards the walls (dashed red 
line). For comparison, the Maxwell Boltzmann distribution 
expected in cm scale vapor cells is shown (dashed black line). 


where / is a narrowing factor for the slow atoms. The 
full distribution is plotted as a black line on Fig. |^b). 
The distribution is then convolved with the collisional- 
broadened line shape. 

To test if this bimodal distribution is appropriate, 
Monte-Carlo (MC) simulation of atomic trajectories was 
performed and the modified position and velocity distri¬ 
butions were included in the line shape model, as well 


Gsiow(i') = exp{-v'^ / (fuf) 
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as a 1/r^ atom surface potential that induces a force 
F = —dUydw/^z across the cell. In this simulation we 
assume that the atom-light interaction region is a box 
with dimensions w xw x L, where w = 40 //m is the 1/e^ 
diameter of the excitation laser beam used in the exper¬ 
iment and L is the length of the nanocell. The atoms 
are randomly placed uniformly across a region which is 4 
times larger than the interaction region {Aw x Aw x L ), 
to account for atoms transiting into the interaction region 
during the simulation. The atoms are given random ve¬ 
locities according to Maxwell-Boltzmann statistics, and 
their motion is simulated over a time (1 //s) which is large 
compared to the lifetime of the excited state. We do not 
include any other decay mechanisms. We neglect desorp¬ 
tion events, so atoms that hit the walls or move out of 
the interaction region are lost. This is justified since des¬ 
orption processes are expected to happen on time scales 
much longer than the interaction time [45] . 

The results of the simulation are shown in Fig.[^ where 
dashed red (blue) lines show the results of the model 
without (with) AS forces present. Panel (a) shows the 
position-dependent interaction times of atoms inside the 
cell. Panel (b) compares the MC simulated (blue line) 
velocity distribution of the atoms to the bimodal dis¬ 
tribution used in the model (solid black line) and the 
Gaussian distribution expected for a Doppler-broadened 
medium (dotted black line). The bimodal distribution is 
a reasonable approximation. 

In panel (d) of Fig. we account for the AS interaction 
by calculating the position-dependent AS shift Aas{t) 
across the entire cell length L: 


Aas 





(L-r)^ 


(9) 


where Ca is a coupling coefficient describing the strength 
of the AS interaction, and a is an exponent describing the 
power law of the AS interaction. The two terms in Eq. 
simply sum the AS interaction from both cell walls. Both 
Ca and a are used as floating fit parameters to determine 
the form of the AS power law. Aas is then weighted by 
the position-dependent interaction times, shown in blue 
in Fig. The result is an asymmetric line shape rep¬ 
resenting shifts across the entire cell weighted by inter¬ 
action times. This line shape is then convolved with the 
Dicke narrowed line shape to give the final result of the 
fitting function. Note that the AS forces slightly modify 
the position and velocity distribution as shown by the 
MC simulation results in Fig. (red dashed lines), but 
this is only significant for distances less than 10 nm. 

Additional broadening effects that are not currently 
modeled are accounted for by adding an excess broaden¬ 
ing parameter Fex into the fit. We suspect these addi¬ 
tional broadenings may be caused by radiation trapping 
within the nanocell, a future avenue of investigation. 

To summarize; the model takes the natural line shape, 
includes the well-known self broadening, and a fitted 
additional broadening, Fex- We then convolve the line 
shape with a phenomenological bimodal velocity distri- 
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FIG. 7. An accumulation of many spectra for different cell 
lengths. Here we highlight the length dependence of the ab¬ 
sorption spectra on the ®^Rb bSi/ 2 Fg = 2 line to the ®^Rb 
5 P 3/2 Fe = 1,2,3 excited state for T = 150°C. The expected 
AS shift for rmax = L/2 is plotted as a white dashed line, using 

Aas = —This data set clearly demon¬ 
strates the length dependence of the red tail caused by the 
AS interaction. However in Rb, the red tail and shifts are not 
appreciable until cell lengths shorter than 80 nm, i.e. maxi¬ 
mum AS separations of 40 nm. Hence, fits performed will not 
be sensitive to the AS interaction outside of this length range. 


bution to account for Dicke narrowing, with a fitted pa¬ 
rameter a representing the strength of the contribution 
of traditional Maxwellian velocity classes. The last step 
is to convolve with a line shape generated using the AS 
interaction, fitting for the AS coefficient for the excited 
state transition Ca and the exponent of the in¬ 

teraction a. 


IV. RESULTS AND DISCUSSION 
A. Rubidium D2 line 

Unlike other experiments that map out the AS inter¬ 
action using peak shift measurements, we choose to in¬ 
vestigate the AS interaction by fitting a full line shape 
which is most sensitive to smaller AS distances. Spectra 
are taken over a length range of L = 40 — 250 nm at 
temperatures T = 125, 150, 175 and 200 °C. The rela¬ 
tive ease of acquisition of spectra also allows us to take 
many readings over long integration times, increasing the 
precision of the experiment. 

A color map showing the len^h dependence of the AS 
induced shifts is shown in Fig. Data are on the ^^Rb 
55 'i/ 2F^ = 2 line to the ^^Rb 5 P 3/2 F^ = 1,2,3 line of the 
Rb D2 resonance taken at T = 150°G. Each horizontal 
slice shows an experimental spectrum like those in Eigs. 

(c) and (d), i.e. with background counts subtracted, 
then normalized such that 1 on the color scale is the 
peak of the ^^Rb 55 'i/ 2 F’^ = 2 to the ^^Rb — I 5 









250 



Cost Value 



FIG. 8. Analysis of simulated annealing (SA) fits for a 
cell of length 80 nm, temperature 125°C. (a) Cost values of 
attempted Ca coefficients tried by the simulated annealing 
fitting algorithm, with the final fitted value plotted as a red 
dashed line. There is no clear best value of Ca over the range 
explored, (b) Cost values of exponent a explored by the SA 
fitting algorithm. Here, there is a clear region between a = 
2.9 and 3.1 where the fit is best. The fitted value found by 
the SA fitting algorithm is shown as a red dashed line, (c) 
Parameter space showing cost value across an array of Ca and 
a with the rest of the fit parameters kept constant. Panels (d) 
and (e) show the cost values along the fit parameters found by 
the SA fitting algorithm, highlighting that the region where 
a is the best fit parameter is very narrow, whereas there is 
smaller relative variation in the cost value for Ca- 


2, 3 excited state transitions. The figure demonstrates 
a clear length dependence of the asymmetric red tail, as 
expected for the AS interaction. The expected AS shift 
for atoms in the center of the cell is plotted as a white 
dashed line. The actual peak shift does not follow exactly 
the expected 1/r^ form because the asymmetric red tail 
skews the peak location. Hence, line-shape fitting should 
be more effective than peak detection for determining the 
functional form of the AS potential. 

At long lengths (L > A/4), shifts due to the AS inter¬ 
action on the ground state transitions are very small - on 
the order of a few MHz. The impact on the line shape 
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FIG. 9. Length dependence of goodness of fit parameter 
T [37], purple points in panel (a), and the fitted AS exponent 
a, red points in panel (b). The results from longer cell lengths 
generally have poorer fits than shorter cell lengths where the 
AS interaction dominates (L < 80 nm), highlighted as a gray 
area. The final value of a = 3.02 ± 0.06 is plotted as a dashed 
black line, with the error bars plotted as beige dashed lines. 


is minimal until shorter cell lengths. For the Rb D2 line 
the asymmetric tail is not appreciable until cell lengths 
less than 80 nm as evidenced in Fig. The asymmetric 
red shift characteristic of the AS interaction is not fully 
evident until the cell length is below 80 nm with maxi¬ 
mum AS separation rmax = 40 nm. This means our fits 
are most effective in the length region r^iax = 10-40 nm, 
where there is a sufficiently large red tail to fit to. 

After background-count subtraction and normaliza¬ 
tion, spectra are fitted using the model outlined in the 
previous section. Examples of fit results on the Rb D2 
line for a temperature T = 125°C are shown in Figs. |^(c) 
and (d). They show experimental spectra (black points) 
processed from the raw data shown in Figs. I^a) and 
(b), fits performed for cell lengths of 60 nm (red line) 
and 250 nm (blue line). 

Upon visual inspection, both spectra have excellent 
fits. The fit for L = 250 nm in panel (b) has an T of 
1.7 indicating an excellent fit, although with no AS in¬ 
duced shifts apparent because of the long cell length. The 
fit for L = 60 nm in panel (a) with an T of 0.1 indicates 
a possible overestimation of the errors involved in the fit- 
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ting procedure m However, as discussed earlier, this is 
caused by the large background counts for such a short 
cell length. 

To extract an exponent a, fit results from all spec¬ 
tra are taken, and an average weighted by the inverse of 
the T of the fit is calculated. The same treatment per¬ 
formed on fitted values of Ca to determine the 

coupling coefficient. Ultimately we have determined the 
best values of a and Ca to use in future fitting 

of absorption spectra. Using this method, we extracted a 
value ofa = 3.02±0.06, a value in good agreement with 
the theoretical expectation of a = 3.16 for the ground 
state interaction, discussed in section in We also find a 
value of 1 4^0.3 kHz /im^, not in agree¬ 
ment with theoretical calculations of =2.1 

kHz yum^ [46]. Deviation of the fitted coefficient from 
theoretical values will be discussed in the next section. 

To test the robustness of our fits, we refitted spectra 
from our data set using a simulated annealing (SA) fitting 
algorithm [52] instead of the Marquardt-Levenberg (ML) 
method. The SA fitting algorithm was adapted from the 
code in [49] to use our model of the absorption spectra 
from section HI To identify the best fit the SA algorithm 
tries fit parameters and calculates a cost value, calculated 
as the absolute value of the sum of the theoretical fit 
subtracted from the dataset. The parameters are varied 
such that a global optimal fit is found as opposed to a 
local optimal fit which can returned by a ML fit. The 
SA algorithm was further adapted to return attempted 
fit parameters along with the cost value returned for that 
attempt. 


Fig.j^shows the resulting cost values for 
panel (a), and a, panel (b). Final fitted values from 

the SA fit, ^ i j^Hz /im^ and a = 2.99 

are plotted as dashed red lines on their respective plots. 
Panel (a) shows that there is very little variation in the 

lowest cost value for (7a over the range of values 

explored by the SA fitting algorithm, indicating a larger 
error bar. Conversely, there is a clear region of lowest 
cost for attempted values of a, centred around the fitted 
value of Of = 2.99, in agreement with previous analysis. 

To explore the possible best fit values, we generated 
a theoretical parameter space. We kept all parameters 
aside from a and (7a constant, then calculated 

the cost values across an array of Ca and a. 

The resulting parameter space is shown in Fig.|^(c). We 
also show the variation in the cost value along the axes 
defined by the fitted SA parameters (purple cross), for 
a in Fig.j^ (d) and for Fig. (e). The 

selections are highlighted on panel (c) as blue dashed 
lines. Again, the region of lowest cost value for a is 
much narrower than for Ca showing that this 

method of fitting spectra is more effective in identifying 
the exponent, a, of the AS interaction. However, the 
larger uncertainty in best cost value for Ca in 

the parameter space shows that this method is less suited 
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FIG. 10. Comparison of the red (a) and blue (b) wings of 
the Cs 6SiFg = 3 to 6PiFe = 4 transition. Zero detuning 
is the resonance frequency of the line. There is a remarkable 
difference between the two wings, showing what a profound 
impact the AS interaction has on the line shape. The data 
(black points) are compared to the red wing as fitted by our 
model (dashed red line, mirrored in panel (b) for comparative 
purposes) and a Lorentzian fit to the blue wing (blue line). 
A Lorentzian fit is sufficient to model the blue wing, because 
Doppler broadening is suppressed by Dicke narrowing of the 
resonance. 


to find an accurate value of 

The length dependence of the T (blue points) of each 
ML fit with the corresponding fitted values of a (red 
points) are plotted in Figs. |^(a) and (b) respectively. In 
panel (b), the experimentally determined value is high¬ 
lighted with a black dashed line, with the standard de¬ 
viation plotted as gray dashed lines. Comparing with 
Fig. i (c), where there is a range of a where the cost 
value is minimal, between a = 3.2 and 2.9, we can see 
that the fitted a from the ML fitting method all lie within 
this length range, showing that the ML fit explored the 
whole region of lowest cost values of a, and still returns fi¬ 
nal fit parameters clustered around our final quoted value 
of a. In panel (a) the T of the fit is generally larger for 
longer cell lengths, where the AS shift is not appreciable. 
Fitted values of a and (7^‘^D2^5P3/2 these 

longer lengths have very little bearing on the nature of 
the AS interaction, and a weighted average using T will 
decrease their impact on the final measured value. For 
shorter cell lengths L < 80 nm, the AS interactions have 
much stronger effects on the line shape also coinciding 
with smaller values of T, with most fitted values of a 
in good agreement with the final average value and its 
associated errors. 


B. Cesium D1 line 

To corroborate the results found in Rb, and to test 
the model with other alkali-metal atoms, the experiment 
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was performed in a Cs NC on the Cs D 1 894 nm line, 
where large hyperfine splittings (9 GHz for the Cs ground 
states and 1.2 GHz on the excited states). Each hyperfine 
transition is separated such that it can be considered in 
isolation from all others, even in the presence of a large 
AS interaction. This allows the reliable measurement of 
shifts with a precision of a few MHz. 

Whilst the ability to address individual transitions 
gives exceptionally high resolution, the technical chal¬ 
lenges of the Cs D 1 line are greater than the Rb D 2 
line, primarily due to the transition wavelength. Firstly, 
the quantum efficiency of Si, used for detection in the 
SPCMs, is much lower at 894 nm than at 780 nm and 
hence the detectors are not as sensitive. Secondly, high 
transmission, narrow bandpass filters are currently not 
available at 894 nm meaning relatively more thermal 
background photons (of which there are more in abso¬ 
lute number compared to 780 nm) reach the detector. 
As a result of these two factors, we limit our investiga¬ 
tion to nanocell thicknesses greater than 100 nm for the 
Cs D 1 line in order to generate spectra with sufficiently 
high signal-to-background ratio. 

Using the same methods as those to extract param¬ 
eters for Rb, we find a value of g ^ 

0.1) kHz /im^, a value that, similar to the Rb coeffi¬ 
cient, is in disagreement with theoretical expectations of 

In Fig. we present a comparison of the red and 
blue wings of AS shifted transitions, the Cs 65'i/2^^ = 3 
to 6 Pi/ 2 ^e = 3, panel (a), and = 4, panel (b). The 
blue line on both plots is a Lorentzian fit, i.e. one that 
excludes the AS interaction. The red line on panel (a) 
is a fit including the AS interaction, and is mirrored in 
panel (b), highlighting the asymmetry of the AS inter¬ 
action. These spectra show that broadening effects are 
exclusively on the red wing. 

V. CONCLUSION 

We measure and analyze the resonance line shapes of 
atoms confined between sapphire plates with separations 
between 30 and 250 nm. The lineshapes are asymmetric 
due to the AS interaction. 

We fit the absorption line shapes in this length range, 
with an AS interaction of the general form —Cajr^. We 
found the best fitted values by fitting absorption spectra 
taken at multiple temperatures over the aforementioned 


length ranges, a method that greatly enhances the preci¬ 
sion of our experiment. The resulting fitted parameters 
were averaged with a weighting according to the reduced 
chi squared (denoted as T) of the fit. We find an averaged 
value of a = 3.02 ± 0.06, a value in good agreement with 
theoretical expectation for the atom surface interaction 
inside the near field, and confirmed by other experiments 

[54]. 

We also found a fitted value of = 14 ^ 

0.1 kHz /im^, a value not in agreement with theoreti¬ 
cal expectations of 2.1 kHz /im^. We postulate that 
the discrepancy between experiment and theory may 
be caused by rubidium adsorbed on the interior sur¬ 
face of the cell, interfering with the surface permittiv¬ 
ity and therefore the reflection coefficient used to calcu¬ 
late [^| 44 ], xhe effect of metallic layers 

on the surface affects only the (.QgfXcient, 

leaving the exponent of unchanged [44] . This leads us 
to conclude that vapour cells may not be an appropri¬ 
ate medium for the accurate measurement of without 
rigorous theoretical treatment to account for the effects 
of surface atoms, currently beyond the scope of this ex¬ 
periment. There may be more complex interactions oc¬ 
curring inside the cell [55], such as additional motional 
effects between photon absorption and fluorescence, mul¬ 
tiple reflections and surface charges. However, we have 
demonstrated that a simple model is sufficient to describe 
the power law of 1 /r^ for spectra over the length scales 
investigated. 

In future we intend to utilize the methods developed 
herein to aid future investigations on the AS interaction, 
including the search for bound states close to the sur¬ 
face [56] and investigate higher-excited D states [17]. To 
further characterize our system, we also hope to look at 
the effects of radiation trapping |57]. Finally, we will use 
the model developed to aid in the fabrication and testing 
of nanocells being developed in Durham. 
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